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ABSTRACT: We recalculate the two-loop beta functions for three gauge couplings taking 
into account all low energy threshold corrections in split supersymmetry (split-SUSY) which 
assumes a very high scalar mass scale Mg. We find that, in split-SUSY with gaugino mass 
unification assumption and a large Ms, the gauge coupling unification requires a lower 
bound on gaugino mass. Combined with the constraints from the dark matter relic density 
and direct detection limits, we find that split-SUSY is very restricted and for dark matter 
mass below 1 TeV the allowed parameter space can be fully covered by XENON-1T(2017). 
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1. Introduction 


It is well known that both the ATLAS and CMS collaborations have established the exis- 
tence of a 125 GeV Standard Model (SM)-like Higgs boson zl So far the LHC Higgs 
data (with large uncertainties) agree well with the SM predictions. Still, such a newly dis- 
covered Higgs boson (especially its enhanced diphoton signal rate reported by ATLAS) has 
been interpreted in various new physics frameworks, among which a particular interesting 


scenario is low energy supersymmetry [B]. 


Supersymmetry (SUSY) is interesting in many aspects. A very interesting observation 
is that the observed Higgs boson mass of 125 GeV falls within the narrow window 115 — 135 
GeV predicted by the Minimal Supersymmetric Standard Model (MSSM). Besides, the 
unification of gauge couplings | H. which cannot be achieved in the SM, can be successfully 
realized by introducing supersymmetric particles with proper quantum numbers. The 
observed cosmic dark matter, which has no interpretation in the SM, can be perfectly 
explained in SUSY. 


Although SUSY is appealing, no signals of SUSY have been found at the LHC, which 
implies that squarks and gluinos should beyond the 1 TeV range. In fact, the LHC data set 
a limit, M mg > 1.5 TeV for mg ~ mg and mg Z 1 TeV for mg > mg within the popular 
CMSSM model. On the other hand, radiative electroweak symmetry breaking conditions 
to give a 125 GeV Higgs requires an electroweak fine-tuning (EWFT). Such a fine-tuning 
may indicate that we should not expect SUSY to provide naturalness. Actually, from the 
viewpoint of quantum field theory, the naturalness problem of the Higgs mass appears 
to be quite similar to the cosmological constant problem, since both of them are related 
to ultraviolet power divergences. Maybe we can apply the naturalness criterion of the 
cosmological constant to SUSY. Split supersymmetry (split-SUSY), proposed in B. p, [L9]. 
gives up naturalness while keeps the other two main virtues: the gauge coupling unification 
and viable dark matter candidates. This split-SUSY scenario assumes a very high scalar 
mass scale Mg and at low energy the supersymmetric particles are only the gauginos and 
higgsinos as well as a fine-tuned Higgs boson. With very heavy sfermions this scenario can 
obviously avoid the flavor problem. 


Given the significant progress of the LHC experiment and dark matter detections 
fill, [12] ER we in this work check the dark matter and gauge coupling unification in split- 
SUSY. In fact, as shown in (4, (15, Ed], the previous dark matter data can already set 
some constraints on the parameter space of split-SUSY. The gauge coupling unification 
in split-SUSY had been checked at two loop level in a special case assuming Mı = My = 
Ma = m and also in complete two loop level in fig}. We recalculate the two-loop 
beta functions for three gauge couplings at two loop level taking into account all threshold 
corrections to check the status of split SUSY after higgs discovery, in particular the gauge 
coupling unification constraints on dark matter phenomenology. 

This paper is organized as follows. In Sec. bl we study the gauge coupling unification 
in split-SUSY. In Sec. [3| we examine the constraints of dark matter relic density and direct 
detections on split-SUSY. Sec. pl contains our conclusions. 


2. Constraints of Split SUSY From Gauge Coupling Unification 


We firstly brief review the split supersymmetry scenario and explain our conventions. More 
details can be found in B. gj. The Lagrangian of Split supersymmetry is given by 
2 2 d 7 
L= mH H-I (HH) — |hsajueH" + Mh qidH + hte H 
Mo ~oz M, zz zp > 
+ 749^ + 22074 + z BB + pH eH, 
Hi ~ aya J Df HTe ~ aya ^! DI] 
De (uo W + 9,B) Ay + Um (-Àuc W° + 34B) Hy+h.c.} , (2.1) 
with e = ic» and the higgsino components Ĥa, the gluino g, the Wino W, the Bino B as 
well as all the standard model particles with one Higgs doublet H. The standard model 
higgs doublet is the linear combination of two higgs doublets H = — cos peH% + sin SA, 
which are fine-tuned to have small mass. The definition of scalar quartic coupling A and 
the yukawa couplings hoe will be given shortly. The parameter u arises from the u-term 
of the supersymmetric standard model and acts as the higgsino mass parameter. 

'The squarks, sleptons, charged as well as the pseudoscalar Higgs from the supersym- 
metric standard model in split SUSY scenario are assumed to be heavy (so that they will 
not cause a problem in SUSY flavor problems etc) and their masses are assumed to be 
degenerated at mass scale Mg. The coupling constants appeared in previous Lagrangian 


at the scale Mg are obtained by matching them with the interaction terms of the super- 
symmetric Higgs doublets H, and Hg 


g9? t.a Ta 2 g’? t 1 2 
Lowy = -E (Hio H, + Hlc Ha) -5 (nin, -H ino) 


ATH, tig; = Ad; Hj ed;qj — Af; HI ceil; 


HS (gorie + of B) Ëu- 8 (got — gB) lath (2.2) 
——= | go g u — — [90 —g athe... . 


Because one Higgs doublet can be fine-tuned to be small, the new coupling constants at 
the scale Mg can be obtained by replacing H, — sin 8H and H4 — cos GeH* into (B.9) 
with: 


[g^ (Ms) + g^ (Ms)) 


A(Mg) = RENI MESES cos? 2f, (2.3) 
ha (Ms) = X (Ms) sin B, h° (Ms) = Ap (Ms) cos B, (2.4) 
gu(Ms) = g(Ms) sin B, ga(Ms) = g(Ms) cos B, (2.5) 
Gu(Ms) = g'(Ms)sin B, ğa (Ms) = g' (Ms) cos B. (2.6) 


We should note that such tree level relation will hold in higher order only if DR (Di- 
mensional Reduction) renormalization scheme is used. Supersymmetry ensures that the 
gaugino coupling g within v28d! (t^)! (à^) is equal to the gauge couping g. Due to the 
fact that MS is not supersymmetry preserving, the relation g = g is spoiled in this scheme. 
The relation (p.3) will be modified to act as the input of RGE running (see appendix). 


Let us take a look at the free parameters in split-SUSY. It is well known that for the 
ratios of gaugino masses and gauge couplings we have 


d (MN. 
— (=) ~0 (2.7) 


and thus the ratios are RGE-invariant (up to one-loop level). This leads to a mass relation 
given by 


Mi M» 


M3 My 
È R GF 


(2.8) 


with universal gaugino mass at the GUT scale. This gaugino mass relation can naturally 
appear in the ordinary SUSY-SU(5) GUT models (it can be spoiled by the introduction 
of certain higher dimensional representation Higgs fields, e.g., the 75, 200 dimensional 
Higgs fields d, Bi). The two-loop corrections to the mass ratios M;/g? are subdominant 
and make negligible contributions to two-loop RGE running of gauge couplings. So in our 
following analysis we adopt this gaugino mass relation. With this mass relation, the low 
energy SUSY mass parameters in split-SUSY can be reduced to: Ma, u and Mg. The 
parameter tan 8 is chosen by random scan so as to give the 125 GeV higgs in the next 
section. It was chosen as a free parameter in this section. To avoid the SUSY flavor 
problem, split-SUSY assumes Ms >> (M3, u) and the value of Ms is typically chosen to 
be higher than 100 TeV. We should note that the gaugino mass relation will no longer be 
valid below Mg due to the split nature of the split supersymmetry spectrum. However, 
various constraints, especially the 125 GeV higgs discovery by LHC, exclude the high Mg 
scenario and favor scalar superpartners in the region Mg ~ 104 — 108GeV [I8]. So it can 
be reasonable to keep the approximate ratio of the gaugino mass relations. 

Preserving gauge coupling unification is one of the two motivations of split-SUSY 
which, on the other side, is a highly non-trivial constraint on split-SUSY. In general, 
the successful gauge coupling unification at one-loop level taking into account threshold 
corrections disfavors a large Ms due to the prediction of a relatively lower a,(Mz) than the 
experimental value. In |8| it is argued that the two-loop renormalization group equation 
(RGE) running can alleviate this difficulty by pushing up the predicted a3(Mz) to around 
0.130 and thus can push up Ms to a large value. So the inclusion of two-loop RGE runnings 
for gauge couplings are necessary in order to achieve the gauge coupling unification in split- 
SUSY. 

In this work we use the method in [2] R3] to calculate the two-loop beta functions for 
three gauge couplings in split-SUSY, taking into account the threshold corrections. The re- 
sults of p, which assuming Mı = Mə = Ms = y, is a special case of our general results (we 
checked that in this special case both results are in agreement). To study the RGE running 
for gauge couplings, we also calculated the one-loop beta functions for Yukawa couplings 
and gaugino couplings with threshold corrections. There are in total four different scenarios 
depending on the relative size of the gaugino masses and u. The full analytic expression 
for the beta function in these scenarios can be seen in the appendix. Although the pro- 
ton decay problem in the split susy scenario will ameliorated, natural doublet-triplet(D-T) 


splitting may still need certain mechanism. Incorporating various D-T splitting mechanism 
can lead to uncertainties in the GUT theory field contents and consequently new matter 
threshold uncertainties. So in our study on gauge coupling unification, we neglect possible 
GUT scale threshold corrections and possible new gauge kinetic terms from Planck-scale 
suppressed non-renormalizable operators involving various high representation higgs fields 
of GUT gauge group. It is well known that the two loop RGE running for gauge couplings 
are scheme independent, so we use the M S couplings in our studying of the gauge coupling 
unification. 

With the two-loop RGE running of gauge couplings, we can study the gauge coupling 
unification requirement for the three free mass parameters in split-SUSY. To make our 
calculation reliable, the GUT scale must be significantly lower than the Planck scale so 
that the gravitational effects can be neglected. On the other hand, the GUT scale can not 
be very low; otherwise it will lead to fast proton decay. 

Note that in ordinary SUSY-GUT, the dominant proton decay comes from the dimension- 
5 operators involving the triplet Higgs and gaugino loops (these dimension-5 operators in- 
duce the decay p + K * --7, whose experimental bound is 7;, , c; > 3.3x 10?? years4, P3). 
Since this decay also involves sfermions in the loops, it is much suppressed in split-SUSY 
due to very heavy sfermions. In fact, as noted in Ed, the contribution from the model- 
dependent dimension-5 operator which is suppressed by M4 is subdominant to dimension-6 
operators if the amplitude is suppressed by two light quark/lepton masses. In Split Su- 
persymmetry, the heavy squarks can provide adequate suppression and the suppression of 
light fermion masses can even be unnecessary. 

So for proton decay, we only consider the decay mode p — et + 7° induced by the 
heavy X, Y gauge bosons of SU(5) with mass Maur through the dimension-6 operators 
(via gauge boson exchange) [p]: 


M * / 1/35 \? /0.015 GeV? / 5 \? 

0.4 GUT 34 
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with Ar, the operator renormalization factors and ay the hadronic matix element. The 
lattice result [dj gives ay = 0.015 GeV? 
Combining with the experimental bound given by[P4. 


T(p > et + 7°) > 1.0 x 1094 years, (2.9) 


we can find the lower limit for the GUT scale. Taking into account the upper limit (Planck 
scale) and choosing the central value of Ar, — 5 in equation Œ), the GUT scale should 
lie in the range 


1.0 x 10!°GeV > Maur > V35acur (6.9 x 107?) GeV . (2.10) 


In our numerical study, we require that successful grand unification should satisfy this 
constraint on the GUT scale. 

The following setting is used in our numerical studies: We use the central value of 
91,92 and 3c range of g3 as the input at the electroweak scale. Other couplings at the 
electroweak scale, for example, the top yukawa h etc, are extracted from the standard 


model inputs taking into account the threshold corrections. Relevant details can be seen 
in the appendix. We also use their central values in our numerical studies. 

Gauge couplings unification requires that the three gauge couplings meet at the same 
point with g1( Maur) = g2(Mcur) = 93(Maur) and the GUT scale satisfied the equation 
(2.10). However, in numerical studies, it is not possible to obtain exact equality which 
differs dramatically from the approach of the one-loop case. Because of the decoupled 
nature of the one-loop gauge couplings running, the unification scale is determined by the 
intersection of g1, go and one can extrapolate back to predict g3 at the elctroweak scale. In 
case of the two loop results, the two-loop RGE running of gauge couplings which amount 
to numerically solve a series of coupled differential equations are obtained from the values 
at electroweak scale and evolve step by step to GUT scale. We thus use the criteria that 
the gauge couplings unification is satisfied when the three couplings differ within the range 
0.005 (less than 1% error). 

The RGE running of the three gauge couplings for some benchmark points in the 
parameter space is displayed in Fig fl} where we fix Ms = 100 TeV, u = 500 GeV, tan 8 = 
10 and vary Mz from 200 GeV to 3.33 TeV. To illustrate if the three gauge couplings can 
really merge at a high scale, we only show the running region of E > 10!4 GeV in this 
figure. In fact, we found that the two-loop RGEs change go coupling more sizably than 
gı and gs. We can see from this figure that gauge coupling unification prefers a relatively 
large gaugino mass. 

With a random scan over the parameter space (0 < Mo, < Ms < 10l?GeV) for 
1 < tan f < 50 under the gauge coupling unification requirement, we obtain the results 
shown in Fig. p) The sharp edge within the figures corresponds to the constraints Ms > M3 
in the split SUSY. From the left panel we can find an upper bound for Ms, which is about 
10° GeV (since split-SUSY requires Mg > Msg,, we can also obtain an upper bound on M» 
correspondingly). From the right panel we can find upper limits for u and M2, which are 
around 100 TeV, independent of the Mg value. 
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Figure 1: The RGE running of the three gauge couplings (we only show the region of E > 1014 
GeV). The dashed lines (green) denote the one-loop results while the solid lines (red) denote the 
two-loop results. 
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Figure 2: The scatter plots of the parameter space with the gauge coupling unification requirement. 
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Figure 3: Same as Fig.2, but showing M2 versus Ms for fixed p. 


We also scan the parameter space of (M2, Ms) with a fixed value of u and display the 
results in Fig B. We can see that the gauge coupling unification imposes a lower bound on 
Msg, which is 5 TeV for a small u value. It is also interesting to note that a lower bound 
for Mo exists for a large u value. However, when u turns small, the lower bound for M» is 
relaxed. 

Note that on the plane of (M2, Mg) the gauge coupling unification requirement gives 
a region instead of a line. The reason is that some uncertainties are involved in gauge 
coupling unification requirement. The first uncertainty comes from the measured gauge 
couplings at Mz scale and in our calculation we considered the 3e range of as( Mz). The 
second uncertainty is that the merging of three gauge couplings at some GUT scale is 
not ’exact’ numerically (in our analysis we require the difference between any two gauge 


couplings to be smaller than 0.005 while the gauge coupling strength is about 0.68). 

We should give a brief comment on the role of parameter tan B in the gauge coupling 
unification. Naively, tan 8 does not appear explicitly in the two-loop gauge coupling beta 
functions. However, tan f can affect the gauge coupling RGE running by showing itself 
in the yukawa couplings and the gaugino couplings g, g’. Numerical studies indicates that 
the unification is not sensitive to the choice of tan B. The parameter M;, u, which define 
the thresholds of gauginos and higgsino, can also affect the gauge coupling unification by 
changing the value of beta functions. 


3. Dark matter in split-SUSY 


In split-SUSY the lightest neutralino Xo is proposed to be the Weakly Interacting Massive 
Particle (WIMP) dark matter candidate. We now check the dark matter issue in split- 
SUSY, using the latest relic density data from Planck and the direct detection limits from 
XENONI100,LUX as well as the future Xeon1T. 

We use the package DarkSUSY PT] to scan the parameter space of split-SUSY in the 
ranges: 


1< tan <50, 0 < (Mo, u) < Ms. (3.1) 


In order to use DarkSUSY to calculate the relic density of dark matter in split susy scenario, 
we use the fact that the effects of heavy sfermions and heavy higgs almost entirely decouple 
when Ms = MA > 5TeV [BJ]. So in our numerical study, we single out the points which 
satisfy the GUT constraints (as that in previous section) and then set Ms = M4 = 10TeV 
in DarkSUSY to carry out dark matter related numerical calculations for such survived 
points. 

In our scan we take into account the current dark matter and collider constraints: 


(1) We use the lightest neutralino X? to account for the Planck measured dark matter 
relic density pj; = 0.1199 + 0.0027 (in combination with the WMAP data (4); 


(2) The LEP lower bounds on neutralino and charginos, including the invisible decay 
of Z-boson; For LEP experiments, the most stringent constraints come from the 
chargino mass and the invisible Z-boson decay. We require that mg+ > 103GeV and 
the invisible decay width l'(Z — xoxo) < 1.71 MeV, which is consistent with the 2c 
precision EW measurement result: I?9"75M < 2.0 MeV. 


inv 


(3) The precision electroweak measurements; 


Indirect constraints from electroweak precision observables such as py, sin? 0l ff and 
My or their combinations (oblique parameters S, T, U)|29]. We require the oblique 
parameters to be compatible with the LEP/SLD data at 20 confidence level Bo}. We 
compute these observables with the formula presented in |81[. 


(4) The combined mass range for the Higgs boson: 123GeV « Mj < 127GeV from 
ATLAS and CMS collaborations of LHC. 


In split-SUSY due to large Ms, log(m/mi) >> 1 will spoil the convergence of the 
traditional loop expansion in evaluating the SUSY effects of Higgs boson self-energy. 
So in order to calculate mass of the SM-like Higgs boson, we use the RGE improved 
effective potential|B2]. 'This computation method is employed in the NMSSM Tools 
package[3]. This package can be applied to the MSSM cases by setting A = kx > 0 
so that the MSSM phenomenology is recovered. 


We calculate the spin-independent (SI) dark matter-nucleon scattering rate with the 
relevant parameters chosen as B4 B3, Bd: i = 0.023, fe = 0.032, He = 0.017, 
i = 0.041 and m aes n — 0.020. In our calculation of the scattering rate, we take 
into account all the contributions known so far (including QCD corrections). For fr, we 
take a more reliable value from the recent lattice simulation B7. 
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Figure 4: The scatter plots of the parameter space for u > 0 satisfying constraints (1-4) including 


dark matter relic density. The triangles (red) cannot achieve the gauge coupling unification. 


In figs A and p, we show the scatter plots of the parameter space satisfying constraints 
(1)-(4) with positive u. In the allowed parameter space, some samples cannot achieve the 
gauge coupling unification, which are marked out with red color in these figures. From 
fig.4, we can see that all the parameter space satisfying constraints (1-4) are excluded by 
GUT constraints for Ms 2 200 TeV. 

We see that the current LUX S and XENON100 direct detection limits are quite 
stringent for split-SUSY, which can exclude a large part of the parameter space allowed by 
other constraints including the dark matter relic density. Note that a strip corresponding 
to a dark matter mass range from 1.0 TeV to 1.3 TeV can survive the combined constraints 
of GUT and dark matter direct detection for Ms < 200 TeV. From a careful analysis we 


found that this strip of parameter space gives a higgsino-like dark matter. Outside this 
strip (i.e. for a dark matter mass below 1 TeV), the survived parameter space can be fully 
covered by the future XENON-1T experiment. In fact, the vast majority of such survived 
parameter spaces had already been excluded by LUX. 

For negative u, the survived parameter spaces are shown in fig | and fig [7 Our nu- 
merical calculations show that in most parameter spaces the results are not very sensitive 
to the sign of u. The minus sign scenario can only revive a very small part of parameter 
spaces which otherwise be excluded in positive u scenario. However, unlike the positive u 
scenario, future XENON-1T experiment is necessary to cover all the survived parameter 
spaces with a dark matter mass below 1 TeV. 

So we can conclude that for a dark matter mass below 1 TeV the split-SUSY under 
current experimental constraints and gauge coupling unification requirement can be fully 
covered by the future XENON-1T experiment. 
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Figure 5: Same as Fig.4, but showing the spin-independent cross section of dark matter scattering 
off the nucleon. The curves denote the limits from LUX and XENON100 as well as the future 
XENON-IT sensitivity. 


4. Conclusion 


We calculated the two-loop beta functions for three gauge couplings in split-SUSY taking 
into account all low energy threshold corrections. In split-SUSY scenario with gaugino 
mass unification assumption and a large Ms, we find that the gauge coupling unification 
requires a lower bound on gaugino mass. Combined with the constraints from the dark 
matter relic density and direct detection limits, we found that split-SUSY is very restricted 
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Figure 6: Same as Fig.4 for u < 0. 


Ms = 200 TeV 


“ee 


250 500 750 1000 1250 250 500 750 1000 1250 


Moy (GeV) 


Figure 7: Same as Fig.5 for u < 0. 


and for dark matter mass below 1 TeV the allowed parameter space can be fully covered 
by XENON-1T(2017). 
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Appendix A: Boundary Value of the RGE Running 


We will use the modified minimal subtraction (MS) scheme in our gauge coupling RGE 
running. 

Taking into account certain threshold contributions, the MS couplings can be ex- 
tracted from the standard model input &œs( Mz) = 0.1184 + 0.0007 by 


^2 

g as(Mz) 

dn Mrs = aima (on) a 
2r 3 n (#75) 


Similarly, we have 


Cem (Mz) 


TH (4.2) 
Oem (M. 16 Mt 
1 + Senile) S In (4) 


Qem(Mz) larg = 


with the Standard Model input ag} (Mz) = 127.916 + 0.015. 
The exact form of effective weak mixing angle in the modified minimal subtraction 
M S scheme is rather complex and we use the given by PDG([B9 
2 ĝ? (Mz ) 


= = (0.23116 + 0.00013. 4.3 
7M) + 9 (M) en 


From the top-quark pole mass Mi|pole = 173.5GeV and taking into account the QCD 
threshold corrections, one-loop electroweak corrections as well as two-loop O(ao;) correc- 
tions, the M S input for top-yukawa coupling is given by [40] 


Mi 
hi( M) = 0.93587 + 0.00557 
(Mt) a (at 


173.15 0.00003 Mn 125 
` ` GeV 


(4.4) 


Mz) — 0.1184 
— 0.00041 E 


0.0007 


In converting the pole top quark mass into MS mass, we neglect the subleading possible 
contributions from gaugino corrections in this stage because of undecided gaugino coupling 


914,24; 91u,2u- 
'The bottom and tau yukawa couplings at Mz scale can be similarly extracted from their 
MS or pole mass m, (MS) = 4.18GeV, mr|poie = 1.776GeV followed by RGE running[7] 


2 12/23 2 —3/80 
g3 23 Mz e^ 80 Mz 
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e e 80. /MzNV 7/9 
h.(Mz) = 0.0102 | 1— —~ 1+ ——I — ; 4.5 
(Mz) ( EN Te 9 n (22) i 


Because of the fact that supersymmetry is not preserved in the MS scheme, the bound- 
ary conditions appeared in (p.3) is valid only in DR scheme and will be spoiled in MS 
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scheme. We know that in case of simple group, the MS gauge couplings are related to the 


DR gauge couplings by the relation [i9] 


UNS = IDE h = CC) (4.6) 


The relation (p.3) in MS scheme will be changed into 


gl (Ms) = g'(Ms) sin B | ke (s n 


g 


du(Ms) = o(Ms)sin [1 s (So? - gr). 


ou 


ga( Ms) = g(Ms) sin 8 p+ + — 162 (žr _ 50") | 


24 
~I M — '(M : 1 1 3 2 1 12 4.7 
a(Ms) = g (Ms)sinB |1 + 16:2 | 82 t g? , (4.7) 


at the Ms scale at tree-level. This result agrees with the results in |41|( and also agrees 
with ref. [17] if we use the tree-level expression c? = g?/(g"? + g?) to eliminate g’). 


At one-loop level, the expression changed into 


Gu( Ms) EM 
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with proper normalization g' = 4/3/5gi. Because such boundary conditions are given at 


the Mg scale while other inputs are given at the weak scale Mz, iterative procedure is 


necessary in the numerical studies. 


Appendix B: Two-Loop RGE for Gauge Couplings in Split Supersymme- 


try 


The 2-loop RGE for SU(3)., SU(2);,U(1)y gauge couplings (93, g2, g1, respectively) are 


given by 


d bi 


dE? (4x)? i 


X Byg- XO dTr(hthe) — dw (G+ 82) - da (82 + 82) | , 
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a=u,d,e 


with the U(1)y normalization g? = 3(gy)? and the relevant coefficients in Table [I]BI3H. 
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32 ly, 21 7 4, 21V p 3 
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The one-loop RGE for Yukawa couplings below the Mg scale can be written as 


165 5 1 = h“ ET FORT el Sy deo Ss $ (ntn = 23] 


d 3 
16n he = nd -acto + AT + c4 S, + ch, So + 5 (nta — iia : 


TU = he ES t CRTT + c$,51 + 65,55 + ji à 
(4.9) 
with 
T = Tr(3h h^ + 3h A3 + Ath’), S, = 5 (a + OD], $2 — : (9i + 92) » 


The relevant coefficients in different scenarios can be found in Table BB 
Upon Ms, we recover the MSSM result and the one-loop RGE for yukawa-type inter- 
actions in the superpotential are 


d 
1652 ZA" = A" |—2ct!9? 437 r(Nt AY) + 3X AY + adtad] , 
165 x! =A E + Tr(3M A + ATAS) + AYTAY + Pd i 


165^ 5 y = ye |—2¢F9? 4 Tr(3rtt ad + ret 2) + Baete l 


(4.10) 
with 
88, 4 T738 , 93 
ch = (Sq AEE ASTA 
The gaugino coupling RGE (upon gaugino, higgsino thresholds and below Mg) can be 
written as 
2d. ur - d 9 3 l2 -9 LN ~/2 ~ ~/ x ~ 
167 aie = —3 Juci gi + Iu — 59u9a + ZIuGu + GaGaGu + Gu(T + es, $1 + cs, $3), 
d . - 5. dus bu PU " 
16r’ Ga = —3óaclg? + qi - 5 Gadi + n2 + Gu. Ua + Ga(T + cs, S1 + cs, 3), 
d. "T 3. 5 3.,. ET x 
161^, = —38,619; + 10a + 599d + FIuIu + 3GGaGu + GL (T + es, $1 + cs, 2), 
d. Lm 3. 3. 3. DEN = 
16n* 94 = 398.9) + 79a + iaje + ağa + 3HGuda + G(T + cs, $1 + cs, 594-11) 
with the coefficient 
d 3 11 ~u,d 3 3 
E = (5979): ce = (5579): CS, = CS» =1, (4.12) 


and the boundary value at Mg scale 

Gu(Ms) = 92(Ms)sin 8, ga(Ms) = 92(Ms) cos B, 

Gu(Ms) = gı(Ms)sin b, g;(Ms) = g1(Ms) cos P. (4.13) 
Below Mə, we can decoupling the effect of wino by setting ĝu = gg = 0. Blow Mj, the 


effect of bino can be decoupled by setting g/, = g/, = 0. Below p, these gaugino interactions 
will decouple. 
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Table 1: The coefficients in two-loop gauge coupling RGE with M3 < u < Mg. 
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